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The two-Equation of State (Two-EoS) model is used to describe the hadron-quark phase transition 
in dense- hot matter formed in heavy-ion collisions. The non-linear Walecka model is used to describe 
the hadronic phase. For the quark phase, the Nambu-Jona-Lasinio model coupled to Polyakov-Loop 
fields (PNJL) is used to include both the chiral and (de) confinement dynamics. The phase diagrams 
are derived from the Gibbs conditions and compared with the results obtained in the Hadron-NJL 
model without confinement. As in the Hadron-NJL case a first order transition is observed, but with 
a Critical- End-Point at much higher temperature, consequence of the confinement mechanism that 
reduces the degrees of freedom of the quark matter in proximity of the phase transition. Particular 
attention is devoted to the phase transition in isospin asymmetric matter. Interesting isospin effects 
are found at high baryon density and reduced temperatures, in fact common also to other quark 
models, like MIT-Bag and NJL model. Some possible observation signals are suggested to probe in 
Heavy-Ion Collision (HIC) experiments at intermediate energies. 
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I. INTRODUCTION 



The exploration of the phase diagram of strongly in- 
teracting matter and search for the signals of the phase 
transition from hadronic to quark-gluon phase are very 
important in both theory and experiment. As a funda- 
mental tool, lattice QCD provides us the best framework 
for investigation of non-pcrturbativc phenomena such as 
confinement and quark-gluon plasma formation at finite 
temperature and vanishing (small) chemical potential [1- 
7]. However, lattice QCD suffers the serious problem of 
the fermion determinant with three color at finite /i. Al- 
though several approximation methods have been taken 
to try to evade this problem [8-12], the validity of lat- 
tice simulations at finite chemical potential is still limited 
to the region fiq/T < 1 [13]. The results obtained with 
fiq/T > 1 should be taken with care. 

On the other hand, many phenomenological models 
[14-17], as well as the more microscopic Dyson-Schwinger 
equations (DSEs) approach [18], have been proposed to 
derive a complete description of QCD phase diagram. 
Among these effective models, the Nambu-Jona-Lasinio 
model (NJL) is a predominant one, since it offers a sim- 
ple illustration of chiral symmetry breaking and restora- 
tion, a key feature of QCD [19-25]. Moreover, it provides 
a complicated phase diagram of color superconductivity 
at high density [26-28]. One deficiency of the standard 
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NJL model is that quarks are not confined. Recently, an 
improved version of the NJL model coupled to Polyakov- 
Loop fields (PNJL) has been proposed [29]. The PNJL 
model takes into account both the chiral symmetry and 
(de) confinement effect, giving a good interpretation of 
lattice data at zero chemical potential and finite temper- 
ature. At the same time it is able to make predictions in 
regions that cannot be presently reached in lattice calcu- 
lation [30-36]. 

Most effective models, including the PNJL model, de- 
scribe the hadron-quark-gluon phase transition based on 
quark degrees of freedom. As a matter of fact, at low tem- 
perature and small chemical potential, QCD dynamics 
should be governed by hadrons. Therefore, it is natural 
to describe the strongly interacting matter with hadronic 
degrees of freedom at low T and small fi and quarks at 
high T and large fi. This picture can be easily real- 
ized following a two equation of state (Two-EoS) model, 
where hadronic and quark phases are connected by the 
Gibbs (Maxwell) criteria. Such approach is widely used 
in describing the phase transition in the interior of com- 
pact star in beta-equilibrium (e.g., [37-43] ). Recently, it 
has also been adopted to explore the phase diagram of 
hadron-quark transition at finite temperature and den- 
sity related to Heavy-Ion Colhsion (HIC) [44, 46-49]. 
Moreover, in these studies more attention was paid to 
isospin asymmetric matter, and some observable effects 
have been suggested to be seen in charged meson yield 
ratio and on the onset of quark number scaling of the me- 
son/baryon elliptic flows in Ref. [46, 47]. Such heavy ion 
connection provides us a new orientation to investigate 
the hadron-quark phase transition, and it can stimulate 
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some new relevant researches in this field. 

We have previously studied the hadron-quark phase 
transition in the Two-EoS model by using the MIT-Bag 
model [47] and NJL model [50] to describe quark mat- 
ter, respectively. In particular a kind of Critical-End- 
Point (CEP) of a first order transition has been found 
at about T = 80 McV and n = 900 MeV when NJL 
model is considered for the quark matter. In this pa- 
per, in order to obtain more reliable results and predict 
possible observables in the experiments, an improved cal- 
culation, within the Two-EoS approach, has been per- 
formed. We take the PNJL Lagrangian to describe the 
properties of quark matter, with the interaction between 
quarks and Polyakov-Loop, where both the chiral and 
(de) confinement dynamics are included simultaneously. 
We are not considering here color pairing correlations, 
that are affecting the isospin asymmetry [49], since in 
heavy ion collisions the high density system will be al- 
ways formed at rather large temperatures [46]. 

We obtain the phase diagrams of hadron-quark-gluon 
phase transition in T — pg and T — fig planes. We 
compare the obtained results with those given in [50] 
where the NJL model is used to describe the quark phase. 
The calculation shows that the phase-transition curves 
are greatly modified when both the chiral dynamics and 
(de) confinement effect are considered, in particular in the 
high temperature and low chemical potential region. We 
still see a first order transition but the CEP is now at 
much higher temperature and lower chemical potential. 
In fact the CEP temperature is much closer to the criti- 
cal temperature (for a crossover) given by lattice calcula- 
tion at vanishing chemical potential. Our results seem to 
stress the importance of an extension of lattice calcula- 
tions up to quark chemical potentials around fiq/Tc ~ 1. 

In addition we address the discussion about the non- 
coincidence of chiral and deconfinement phase transition 
at large chemical potential and low temperature, relevant 
to the formation of quarkyonic matter. 

Finally the calculation confirms that the onset den- 
sity of hadron-quark phase transition is much smaller in 
isospin asymmetric than that of symmetric matter, and 



where the antisymmetric tensors of vector mesons are 
given by 

w^!. = - dyUJ^j,, p^^ = d^p^ - dup^. 
The nucleon chemical potential and effective mass in 



therefore it will be more easy to probe the mixed phase 
in experiments. 

The paper is organized as follows. In Section II, we 
describe briefly the Two-EoS approach and give the rele- 
vant formulae of the hadronic non-linear Walecka model 
and the PNJL effective theory. In Section III, we discuss 
the expected effects of the confinement dynamics. The 
quark matter phase transition are presented in Section 
IV for the NJL as well as the PNJL models. Section V is 
devoted to the phase diagrams within the Two-EoS frame 
to the comparison with the results using only the pure 
quark PNJL model to describe both phases. Moreover, 
we present some discussions and conclusions about the 
phase transition, as well as some suggestions for further 
study. Finally, a summary is given in Section VI. 

II. HADRON MATTER, QUARK MATTER 
AND THE MIXED PHASE 

In our Two-EoS approach, the hadron matter and 
quark matter are described by the non-linear Walecka 
model and by the PNJL model, respectively. For the 
mixed phase between pure hadronic and quark matter, 
the two phases are connected each other with the Gibbs 
conditions deduced from thermal, chemical and mechan- 
ical equilibriums. In this section, we will first give a 
short introduction of the nonlinear Walecka model for 
the hadron matter and the PNJL model for quark mat- 
ter, then we construct the mixed phase with the Gibbs 
criteria based on baryon and isospin charge conservations 
during the transition. 

For hadron phase, the non-linear Relativistic Mean 
Field (RMF) approach is used, which provides an ex- 
cellent description of nuclear matter and finite nuclei as 
well as of compressed matter properties probed with high 
energy HIC [44, 46, 47, 51, 52]. The exchanged mesons 
include the isoscalar-scalar meson a and isoscalar-vector 
meson lu [NL force, for isospin symmetric matter), 
isovector-vector meson p and isovector-scalar meson 5, 
{NLp and NLpS forces, for isospin asymmetric matter). 

The effective Lagrangian is written as 
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nuclear medium can be expressed as 

f^i = ^h + <)u>^ + 9pT3iP , (2) 

and 

M* = M - g^a - g^r.^S, (3) 
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where M is the free nucleon mass, t. 



ton and r. 
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-1 for neutron, and /i* 



= 1 for pro- 
is the effec- 
tive chemical potential which reduces to Fermi energy 

at zero temperature. The baryon 



E 



Fi 



M* 



and isospin chemical potentials in the hadron phase are 
defined as 
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The energy density and pressure of nuclear matter at 
finite temperature are derived as 
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where fi{k) and /i(fc) are the fermion and antifermion 
distribution functions for proton and neutron (i ~ p, n): 
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The effective chemical potentials fi* are determined by 
the nucleon densities 



(fik 



(Mk) - Mk)). 



(9) 



With the baryon number density p = Pb = Pp + Pn and 
isospin density = Pp — Pn- The asymmetry parameter 
can be defined as 
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In this study the parameter set NLpS [46] will be used 
to describe the properties of hadron matter. The model 
parameters is determined by calibrating the properties 
of symmetric nuclear matter at zero temperature and 
normal nuclear density. Our parameterizations are also 
tuned to reproduce collective flows and particle produc- 
tion at higher energies, where some hot and dense matter 
is probed, see [52] and refs. therein. 

We take the PNJL model to describe the quark matter. 
In the pure gauge theory, the Polyakov-Loop serves as an 
order parameter for the Z3 symmetry breaking transi- 
tion from low to high temperature, i.e. for the transition 
from a confined to a deconfined phase. In the real world 
quarks are coupled to the Polyakov-Loop, which explic- 
itly breaks the Z3 symmetry. No rigorous order parame- 
ter is established for the deconfinement phase transition. 



However, the Polyakov loop can still be practicable to 
distinguish a confined phase from a deconfined one. 

The Lagrangian density in the three-flavor PNJL 
model is taken as 
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where q denotes the quark fields with three flavors, u, d, 
and s, and three colors; mo = diag(m„, m^, rris) in 
flavor space; G and K are the four-point and six-point 
interacting constants, respectively. The four-point in- 
teraction term in the Lagrangian keeps the S'J7y(3) x 
SUa{^) X Uv{l) X Ua{1) symmetry, while the 't Hooft 
six-point interaction term breaks the C/a(1) symmetry. 

The covariant derivative in the Lagrangian density is 
defined as = 9^ — iAf^. The gluon background field 
Afj^ = S'^Aq is supposed to be homogeneous and static, 

with ^0 = 5-^0 V' '^h^^^ V are SU{S) color generators. 
The effective potential iY($[A], $[A], T) is expressed in 
terms of the traced Polyakov loop $ = {TicMZ-^c and 
its conjugate $ = {TicM) /Nc- The Polyakov loop L is 
a matrix in color space 



L{x) = Vexp 



drAi {x, t) 



(12) 



where /? = 1/T is the inverse of temperature and A4 = 
iAo. 

The Polyakov loop can be expressed in a more intuitive 
physical form as 



$ = exp[-/3Fg(f)] 
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where Fq is the free energy required to add an isolated 
quark to the system. So it wiU go from zero in the con- 
fined phase up to a finite value when deconfinement is 
reached [53]. 

Different effective potentials are adopted in the litera- 
ture [30, 54, 55]. The logarithmic one given in [54] will 
be used in our calculation, which can reproduce well the 
data obtained in lattice calculation. The corresponding 
effective potential reads 



iY($,$,r) 
2^4 



a(T) 
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-6(r)ln 
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gauge sector. Tq is found to be 270 MeV as the criti- 
cal temperature for the deconfinement phase transition 
of the gluon part at zero chemical potential [56]. When 
fermion fields are included, a rescaling of Tq is usually im- 
plemented to obtain consistent result between model cal- 
culation and full lattice simulation which gives a critical 
phase-transition temperature Tc = 173 ±8 MeV [1, 2, 4]. 
In this study we rescale Tq = 210 Mev so as to produce 
Tc ~ 171 MeV for the phase transition temperature at 
zero chemical potential. 

In the mean field approximation, quarks can be seen as 
free quasiparticles with constituent masses Mi, and the 
dynamical quark masses (gap equations) are obtained as 



where 



and 



a{T) = ao + ai ( Y 
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We note that in this version of PNJL the direct cou- 
pling between quark condensates and Polyakov loop is 
only via the covariant derivative in the Lagrangian den- 
sity Eq.(ll). 

The parameters a^, bi are precisely fitted according 
to the lattice result of QCD thermodynamics in pure 
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with i = u,d, s, where 4>i stands for the quark conden- 
sate. The thermodynamic potential of the PNJL model 
at the mean field level is expressed as 
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Here ^{p) = -{if - AU + loifJ-i - «^4)), with fj.^ quark 
chemical potential, is the inverse fermion propagator in 
the background field A4, and the trace has to be taken in 
color, flavor, and Dirac space. After summing over the 
fermion Matsubara frequencies, p'^ ~ iujn = (2n + 1)itT, 
the thermodynamic potential can be written as 
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where Ei = ^Jp"^ + Mf is the energy of quark flavor i. 

We remark some interesting differences with respect to 
the thermodynamical potential derived within the pure 
NJL model, see [24, 50]. Apart the presence of the effec- 
tive potential Z^($, $, T), the Polyakov loop is mostly act- 
ing on the quark-antiquark distribution functions, in the 
direction of a reduction, on the way to conflnement. This 
is largely modifying the quark pressure, as seen in the 
calculations. Moreover, in spite of the minimal coupling 
introduced in the Lagrangian Eq.(ll), only in the covari- 
ant derivative, the quark condensates will be strongly 
affected by the Polyakov loop via the modified g, q dis- 
tribution functions. This will be also clearly observed in 



the comparison of NJL and PNJL phase diagrams. 

The values of cjju, (I'd, it's, ^ and $ are determined by 
minimizing the thermodynamical potential 
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All the thermodynamic quantities relevant to the bulk 
properties of quark matter can be obtained from Q. Es- 
pecially, the pressure and energy density should be zero 
in the vacuum. 

The baron (isospin) density and baryon (isospin )chem- 
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ical potential in quark phase are defined as follows 
1 



= 3(P" Pi^Pu-pd, (21) 



and 



3 1 

Ms = ^(Mu + P-s ^ -^{f^u - f^d)- (22) 

The corresponding asymmetry parameter of quark phase 
is defined as 



p5 Pu + Pd 



(23) 



As an effective model, the (P)NJL model is not renor- 
malizeable, so a cut-off A is implemented in 3-momentum 
space for divergent integrations. We take the model pa- 
rameters: A = 603.2 MeV, GA^ = 1.835, KK^ = 12.36, 
"niu^d = 5.5 and = 140.7 MeV, determined by fitting 
niK and m,, to their experimental values [25] . 
The coefficients in Polyakov effective potential are listed 
in Table I. 



TABLE I: Parameters in Polyakov effective potential given 
in [54] 



ao 


ai 


a2 


as 


3.51 


-2.47 


15.2 


-1.75 



So far we have introduced how to describe the hadronic 
and quark phase by the RMF hadron model and PNJL 
quark model, respectively. The key point in the Two-EoS 
model is to construct the phase transition from hadronic 
to quark matter. As mentioned above, the two phases 
are connected by Gibbs criteria, i.e., the thermal, chemi- 
cal and mechanical equilibrations being required. For the 
Hadron-quark-gluon phase transition relevant to heavy- 
ion collision of duration about 10~^^sec, (10 — 20 fm/c), 
thermal equilibration is only possible for strongly in- 
teracting processes, where baryon number and isospin 
conservations are preserved. So the strange-antistrange 
quark number may be rich, but the net strange quark 
number should be zero before the beginning of hadroniza- 
tion in the expansion stage [57], which can be approxi- 
mately realized by requiring fig = (Hadronization is out 
of the range of this study ) . 

Based on the conservations of baryon number and 
isospin during strong interaction, the Gibbs conditions 
describing the phase transition can be expressed by 



Pb{Pb^Pz^T) = h^{pb,p^,T) 
p"{Pb,P3iT) = fi'i^ipB,P3,T) 
P''{pB,P3,T)^P'i{pB,P3,T), 



(24) 



where pg = {1 ~ x)Pb + XP% and p^ = {1 ~ x)p" + XP^ 
are the total baryon density, isospin density of the mixed 



phase, respectively, and x is the fraction of quark matter. 
The global asymmetry parameter a for the mixed phase 
is 



_P3_ _ (1 - X)p" + XP? 
Pb 



(1 



x)pf + XP% 



(25) 

which is determined by the heavy-ion source formed in 
experiments. 
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FIG. 1: Pressure of quark matter as function of baryon den- 
sity at different temperatures in the NJL model. Isospin sym- 
metric matter. In the shaded area we show also the Hadron 
(NL) curves in the temperature region between 75 and 100 
MeV. 




FIG. 2: Pressure of quark matter as function of baryon den- 
sity at different temperatures in the PNJL model. Isospin 
symmetric matter. In the shaded area we show also the 
Hadron (NL) curves in the temperature region between 150 
and 170 MeV. 
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III. EXPECTED EFFECT OF THE 
CONFINEMENT DYNAMICS 

Before showing detailed phase diagram results within 
the Two-EoS approach it is very instructive to analyze 
the effects of chiral and (dc) confinement dynamics in the 
pure quark sector. In order to understand the physics 
which is behind we will show separately the results in 
the NJL, with the same parameters given before, and the 
PNJL model, for isospin symmetric matter. In Figs. 1, 2 
we plot the pressure of isospin symmetric quark matter 
as function of baryon density for different temperatures 
respectively for the NJL and PNJL models. In this cal- 
culation isospin symmetric matter with IJ. ~ fJ-u — fJ-d 
and /^s = is considered and (f>i stands for the chiral 
condensate oi u,d quarks. 

From the two figures, we can see that the pressure has 
a local maximum and a local minimum at low tempera- 
ture. The local extrema will disappear with the increase 
of temperature. The temperature with the disappearance 
of the two local extrema corresponds to the Critical-End- 
Point CEP of the first order chiral transition, for a more 
detailed discussion please refer to [24, 31]. It is inter- 
esting to note that the critical temperature of the chiral 
transition is rather different in the two cases, around 70 
MeV in the NJL and around 130 MeV in the PNJL, while 
the density region is not much affected. This is due to 
the fact that for a fixed baryon density (or chemical po- 
tential) the NJL presents a much larger pressure for a 
given temperature, as clearly seen from the two figures 
1, 2 [58]. This is a nice indication that when we have 
a coupling to the deconfincment, even if in the minimal 
way included here, the quark pressure at finite tempera- 
tures is reduced since the quarks degrees of freedom start 
to decrease. 

All that will imply important differences at higher tem- 
peratures since above the chiral restoration the quark 
pressure will rapidly increase reaching an end-point in 
the Two-EoS approach where the matching to the hadron 
pressure will not be possible. This will happen in differ- 
ent points of the (T, fi), (T, p) planes for the Hadron-NJL 
[50] and the Hadron-PNJL, and higher temperatures will 
be requested in the PNJL case. In fact this can be also 
clearly seen from Figs. 1 and 2, of the NJL- and PNJL- 
pressures, where we plot also the corresponding curves 
of the hadronic EoS in the end-point regions (shaded 
area). The Gibbs (Maxwell) conditions have no solution 
with decreasing density (chemical potential) and increas- 
ing temperature if we encounter a crossing of the hadron 
and quark curves in the T — ps {T — fj^s) plane, with the 
quark pressure becoming larger than the hadron one. Wc 
see that this is happening for T ~ 75 MeV and p/ po — 1-8 
in the NJL case and for T ~ 170 MeV and p/po — 1.6 in 
the PNJL quark picture. In conclusion, due to the notice- 
able quark pressure difference at finite temperatures, be- 
sides the Critical-End-Points, we expect in general rather 
different phase diagrams given by the Hadron-NJL and 
Hadron-PNJL models. This will be seen in the Section 



V, Figs. 8 and 9. 
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FIG. 3: Chiral condensate (pi (normalized to the vacuum 
value) and Polyakov-Loop $ as functions of temperature for 
various values of the quark chemical potential. Isospin sym- 
metric matter. 



IV. PNJL PHASE DIAGRAM IN THE QUARK 
SECTOR 

In order to better understand the effects of the coupling 
between quark condensates and Polyakov loop and also 
to compare with the Two-EoS results, we discuss here 
also the Phase diagram in the pure quark sector obtained 
from the PNJL model. 

Wc plot in Fig. 3 the temperature evolution of the chi- 
ral condensate </>/ and the Polyakov-Loop $ for various 
values of the quark chemical potential. <& and $ have the 
same values at zero chemical potential and their differ- 
ence is very small at finite chemical potential, hence we 
only present the results of $ in Fig. 3 and later in the 
discussion. 

Firstly, we can see that the chiral condensate and 
Polyakov loop $ vary continuously at ^ = and 200 MeV, 
and there exist sharp decreases (increases) at high tem- 
perature indicating the onset of chiral (deconfincment) 
phase transitions. These characteristics show that the 
corresponding chiral and deconfincment phase transitions 
are crossovers for small chemical potential at high tem- 
perature. At variance, for large chemical potentials, e.g., 
p = 350 MeV, the chiral condensate varies discontinu- 
ously with the temperature, which indicates the presence 
of a first order chiral phase transition, as already seen in 
the pure NJL approach, although at much lower tem- 
perature [24], as discussed in the previous Section. The 
Polyakov loop is always showing a continuous behavior 
indicating that we have only crossover transitions. The 
jump observed for the dash-dotted curve corresponding 
to a. p — 350 MeV chemical potential is just an effect 
of the coupling to the sharp variation of the quark con- 
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densates at the first order chiral transition. Moreover 
this is happening in a region of very small values of the 
$ field at lower temperatures. As a matter of fact such 
discontinuity disappears for the results at /i = 400 MeV, 
i.e. above the chiral transition, see the dashed curves for 
both (j) and $ fields. 
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FIG. 4: Phase diagram of the PNJL model. The correspond- 
ing chiral phase transition for the pure NJL model is also 
shown. Isospin symmetric matter. 



Finally in Fig. 4 we plot the phase diagram of the 
PNJL model inT — /iq plane (always for isospin symmet- 
ric matter). The phase transition curves are obtained 
by requiring d(j)i/dT and d^/dT taking the maximum 
values. For deconfinement phase transition the use of 
the maximum value of d^/dT as the phase-transition 
tracer is a good choice when fj, is not too large. In fact 
we see from Fig. 3 a sharp increase of $ for /i = and 
200 MeV. However, with the increase of chemical poten- 
tial, although we are still able to observe the maximum 
of d^/dT, the width of the maximum increases. The 
peaks of d^/dT arc more smoothed and this will be not 
any more a well defined phase-transition parameter as fi 
is large. Therefore, some authors take $ ~ 1/2 as the 
phase transition parameter [55, 59, 60]. In conclusion 
from Fig. 3 we can see that the chiral phase transition 
is continuous at high temperature and relatively smaller 
chemical potential, while a first order phase transition 
takes place at low temperature and larger chemical po- 
tential. The Critical-End-Point {CEP) of the chiral tran- 
sition appears at (132.2, 296.6) MeV in the T ~ fiq plane, 
in agreement with similar calculations [60]. At variance, 
the deconfinement phase transition is always a continuous 
crossover in the PNJL model, but the peak of 9$ /dT be- 
comes more and more smooth with the increase of baryon 
chemical potential. In addition, at large chemical po- 
tential, a chirally restored but still confined matter, the 
quarkyonic matter, can be realized in the PNJL model. 
All that is reported in Fig. 4 where we plot the full phase 



diagram of the PNJL approach. 

Here we give a short discussion about the coincidence 
of chiral and deconfinement phase transitions as well 
as the presence of quarkyonic matter. The tempera- 
ture dependence of the chiral condensate and of the 
Polyakov loop of Fig. 3 as well the PNJL phase dia- 
gram of Fig. 4 are obtained with the rescaled param- 
eter To — 210 MeV. The coincidence of chiral restora- 
tion and deconfinement takes place at about jjLq = 290 
MeV. If we take Tq = 270 MeV, the approximate coin- 
cidence, with the different phase-transition temperatures 
less than 10 MeV, will move down to /x ~ 0. In any 
case, there is only one cross point of the two phase tran- 
sitions. Up to now, the relation between chiral-symmetry 
restoration and deconfinement phase transition is still an 
open question. It is possible that the coincidence of chi- 
ral and deconfinement phase transition takes place in a 
wider range of chemical potentials. Such coincidence in- 
deed has been recently realized by considering a larger 
coupling [entanglement) between chiral condensate and 
Polyakov loop, with an explicit ^-dependence of the con- 
densate couplig G($) and a chemical potential dependent 
To [59, 60]. 

In the same Fig. 4 we report also the chiral transi- 
tion curve for the pure NJL model (same parameters). 
We note the the coupling between the chiral condensates 
and the Polyakov-Loop fields ($, $) is mostly affecting 
the temperature of the chiral CEP as expected from the 
pressure discussion of the previous Section. 

From Fig. 4 we can see that the deconfinement phase- 
transition temperature is still high at large chemical po- 
tential, and so the region of quarkyonic matter appears 
very wide. On the other hand the signature of a de- 
confinement transition is disappearing at large chemical 
potentials and lower temperatures. Because of the lack 
of lattice QCD data at large real chemical potentials, 
more investigations are needed to study the physics in 
this range. The results in [59] also show that the range 
of quarkyonic matter shrinks when a ^-dependent To 
and /or a larger entanglement between quark condensate 
and Polyakov loop is considered. 

We remark that this (T — /i) zone just represents the 
nuclear mettcr phase diagram region possibly reached in 
the collision of heavy ions at intermediate energies and 
so it is of large interest to perform Two-EoS predictions, 
wich should have a good connection to the more fun- 
damental results of effective quark models. This is the 
subject of the next Section. 



V. HADRON-QUARK PHASE TRANSITIONS 

In the following we will discuss the phase diagrams ob- 
tained in the Two-EoS model, i.e., explicitly considering 
a hadronic EoS with the parameter set of NL for sym- 
metric matter and NLpS for asymmetric matter at low 
density and chemical potential [45, 46]. 

We present firstly the phase transition from hadronic 
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FIG. 5: Phase diagram inT — pg plane in the Two-EoS model 
for symmetric matter. 



to deconfined quark phase in T — pg plane in Fig. 5 for 
symmetric matter and in Fig. 6 for asymmetric matter 
with the global asymmetry parameter a — 0.2. 
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FIG. 6: Phase diagram in T — pg plane in the Two-EoS model 
for asymmetric matter with the global asymmetry parameter 
a = 0.2. X represents the fraction of quark matter. 



For symmetric matter at a fixed T, the first order phase 
transition takes place with the same pressure and pg in 
both phases, but a jump of pg to just as shown in 
Fig. 5. In the mixed phase, the pressures of both phases 
keep unchanged and a = = a9 = for any quark 
fraction x- These features are quite different for the 
mixed phase in isospin asymmetric matter. As already 
noted in [50], where the NJL quark EoS has been used, 
also in the PNJL case we see a clear Isospin Distillation 
effect, i.e., a strong enhancement of the isospin asymme- 



try in the quark component inside the mixed phase, as 
reported in Fig. 7, where the asymmetry parameter in 
the two components are plotted vs. the quark fraction x- 
As a consequence the pressure in the mixed phase keeps 
rising with more rapidly for quark concentrations be- 
low 50 % [50]. 

From Fig. 7 we remark that this isospin enrichment of 
the quark phase is rather robust vs. the increasing tem- 
perature. This is important since color pairing correla- 
tions at low temperatures will decrease symmetry energy 
effects [49]. We have to note that such large isospin dis- 
tillation effect is due to the large difference in the sym- 
metry terms in the two phases, mainly because all the 
used quark effective models do not have explicit isovec- 
tor fields in the interaction [47]. 
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FIG. 7: The behavior of local asymmetric parameters 
and ofi in the mixed phase for several values of temperature. 
Parameter set NLpS is used in the calculation. 



Such behavior of the local asymmetry parameters will 
possibly produce some observational signals in the follow- 
ing hadronization during the expansion. We can expect 
an inverse trend in the emission of neutron rich clusters, 
as well as an enhancement of 7r~/7r"'", /K^ yield ratios 
from the high density n-rich regions which undergo the 
transition. Besides, an enhancement of the production 
of isospin-rich resonances and subsequent decays may be 
found. For more details one can refer to [47, 50]. More- 
over, an evident feature of Fig. 6 is that the onset density 
of hadron-quark phase transition for asymmetric matter 
is much lower than that of the symmetric one, and there- 
fore it will be easier to probe in heavy-ion collision ex- 
periments. 

We plot the T — pg phase diagrams in Fig. 8 for sym- 
metric matter and Fig. 9 for asymmetric matter. Fig. 8 
clearly shows that there is only one phase-transition 
curve in the T — fj,g plane. The phase transition curve 
is independent of the quark fraction x- However, for 
asymmetric matter, the phase transition curve varies for 
different quark fraction x- The phase transition curves 
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in Fig. 9 are obtained with x = and 1, representing the 
beginning and the end of hadron-quark phase transition, 
respectively. 

In Fig. 8 and Fig. 9 we also plot the phase transi- 
tion curves with the Hadron-NJL model. For the NJL 
model with only chiral dynamics, no physical solution ex- 
ists when the temperature is higher than ^ 80 MeV. The 
corresponding temperature is enhanced to about ^ 166 
MeV with the Hadron-PNJL model, which is closer to 
the phase transition (crossover) temperature given by 
full lattice calculation at zero or small chemical poten- 
tial [1. 2, 4]. In this sense the Hadron-PNJL model gives 
significally different results and represents certainly an 
improvement respect to the Hadron-NJL scheme of ref. 
[50]. From Fig. 9 we remark that in both cases the region 
around the Critical — End — Points is not affected by 
isospin asymmetry contributions, which are relevant at 
lower temperatures and larger chemical potentials. 
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FIG. 8: Phase diagram in T — plane for symmetric matter. 



From the detailed discussions of the previous two Sec- 
tions now we nicely understand the large difference be- 
tween Hadron-NJL and Hadron-PNJL phase transitions 
and the important role of the confinement dynamics. 

Finally in Fig. 10 we present together the phase di- 
agrams obtained by the PNJL model and the Hadron- 
PNJL model. We find that the deconfined phase transi- 
tion curve in the PNJL model is close to that obtained in 
the Hadron-PNJL model at high temperature and inter- 
mediate chemical potential. At larger chemical potential, 
the deconfinement phase transition curve in the PNJL 
model has still a high temperature. On the other hand 
from the previous Section we have seen that deconfine- 
ment phase transition order parameter $ cannot describe 
well the phase transition at larger chemical potential and 
lower temperatures. We must rely on the predictions of 
the Two-EoS approach, which in fact nicely show a good 
connection to the results more reliable of the PNJL quark 
model, at high temperature and small or vanishing chcm- 
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FIG. 9: Phase diagram in T — /ig plane for asymmetry matter 
with the global asymmetry parameter a = 0.2. 
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FIG. 10: Phase diagrams of the PNJL model and the Two- 
EoS model (dashed curve). The shaded area is just a guide 
for the eye. 



ical potential. The Two-EoS Hadron-(P)NJL model also 
shows that the phase transition at low temperature takes 
place at much larger chemical potential, consistent with 
the expectation of a more relevant contribution from the 
hadron sector [51]. 

Wc notice that at T = there is no difference between 
the Hadron-NJL and Hadron-PNJL models. This is due 
to the fact that there is no dependence of $ on /xb, there- 
fore it vanishes and the PNJL reduces to the NJL. This 
may casts some doubts on the reliability of the present 
calculations at T = and large fis- However our main 
interest is a region at finite T ( T ~ 50 - 100 MeV) and 
{ fJ'B — 1000 — 2000 MeV) region that can be reached 
by Heavy Ion Collisions at relativistic energies. 

Moreover the results obtained by the Hadron-PNJL 
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model at high T, smah fig and low T, large may be 
improved with the consideration of a stronger entangle- 
ment between chiral condensate and Polyakov loop, and 
a chemical potential dependent Tq [59, 60]. The relevant 
investigation will be performed further study. In 
any case, since we lack of reliable lattice data at large 
chemical potential, in general more theoretical work is 
encouraged. 

VI. SUMMARY 

In this study, the hadron-quark phase transition are in- 
vestigated in the Two-EoS model. The nonlinear Walecka 
model and the PNJL(NJL) model are used to describe 
hadron matter and quark matter, separately. We fol- 
low the Gibbs criteria to construct the mixed phase with 
baryon number and isospin conservations, likely reached 
during the hot and dense phase formed in heavy-ion col- 
lision at intermediate energies. The parameters in both 
models are well fitted to give a good description of the 
properties of nuclear matter (even isospin asymmetric), 
at saturation as well as at higher baryon densities, or lat- 
tice data at high temperature with zero/small chemical 
potential. 

The phase diagrams for both symmetric and asymmet- 
ric matter are explored in both T — p^ and T — fi^ planes. 
In both Hadron- (P)NJL calculations we get a first order 
phase transition with a Critical-End-Point at finite tem- 
perature and chemical potential. In the PNJL case the 
C EP is shifted to larger temperatures and smaller chem- 
ical potential, to the (166, 600) MeV point in the (T, pb) 
plane. This appears a nice indication of a decrease of 
the quark pressure when confinement is accounted for. 
Such result is particularly interesting since the CEP is 
now in the region of fiq/Tc — 1 (where pq is the quark 
chemical potential) and so it could be reached with some 
confidence by lattice-QCD complete calculations. 



Another interesting result is that isospin effects are al- 
most negligible when we approach the CEP. At variance 
the calculation shows that the onset density of asym- 
metric matter is lower than that of symmetric matter. 
Moreover in the mixed phase of asymmetric matter, the 
decreasing of local asymmetry parameter and a*^ 
with the increasing quark fraction x may produce some 
observable signals. In particular we remark the notice- 
able isospin distillation mechanism (isospin enrichment 
of the quark phase) at the beginning of the mixed phase, 
i.e. for low quark fractions, that should show up in the 
hadronization stage during the expansion. We also see 
from Fig. 7 that this effect is still there even at relatively 
large temperatures, certainly present in the high den- 
sity stage of heavy ion collisions at relativistic energies 
[46, 52]. All that support the possibility of an experimen- 
tal observation in the new planned facilities, for example, 
FAIR at GSI-Darmstadt and NICA at JINR-Dubna, with 
realistic asymmetries for stable/unstable beams. Some 
expected possible signals are suggested. 

Because of the lack of lattice data at larger real chem- 
ical potential, we are left with the puzzle between chiral 
symmetry restoration and deconfincment. More investi- 
gations on the chiral dynamics and (de) confinement, as 
well as their entanglement are needed. The improvement 
of the understanding of quark-matter interaction is bene- 
ficial to get more accurate results in the Two-EoS model. 
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